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Abstract 

We study random surfaces which arise as height functions of random perfect match- 
ings (a.k.a. dimer configurations) on an weighted, bipartite, doubly periodic graph G 
embedded in the plane. We derive explicit formulas for the surface tension and lo- 
cal Gibbs measure probabilities of these models. The answers involve a certain plane 
algebraic curve, which is the spectral curve of the Kasteleyn operator of the graph. 
For example, the surface tension is the Legendre dual of the Ronkin function of the 
spectral curve. The amoeba of the spectral curve represents the phase diagram of the 
dimer model. Further, we prove that the spectral curve of a dimer model is always a 
real curve of special type, namely it is a Harnack curve. This implies many qualitative 
and quantitative statement about the behavior of the dimer model, such as existence of 
smooth phases, decay rate of correlations, growth rate of height function fluctuations, 
etc. 
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1 Introduction 

A perfect matching of a graph is a collection of edges with the property that each vertex is 
incident to exactly one of these edges. A graph is bipartite if the vertices can be 2-colored, 
that is, colored black and white so that black vertices are adjacent only to white vertices 
and vice versa. 

Random perfect matchings of a planar graph G — also called dimer configurations — are 
sampled uniformly (or alternatively, with a probability proportional to a product of the 
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corresponding edge weights of G) from the set of all perfect matchings on G. These so-called 
dimer models are the subject of an extensive physics and mathematics literature. (See jH] 
for a survey.) 

Since the set of perfect matchings of G is also in one-to-one correspondence with a class of 
height functions on the faces of G, we may think of random perfect matchings as (discretized) 
random surfaces. One reason for the interest in perfect matchings is that random surfaces 
of this type (and a more general class of random surfaces called solid- on- solid models) are 
popular models for crystal surfaces (e.g. partially dissolved salt crystals) at equilibrium. 
These height functions are most visually compelling when G is a honeycomb lattice. In this 
case, we may represent the vertices of G by triangles in a triangular lattice and edges of G 
by rhombi formed by two adjacent triangles. Dimer configurations correspond to tilings by 
such rhombi; they can be viewed as planar projections of surfaces of the kind seen in Figure 
^ The third coordinate, which can be reconstructed from the dimer configuration uniquely, 
up to an overall additive constant, is the height function. 




Figure 1: On the left is the height function of a random volume-constrained dimer configu- 
ration on the honeycomb lattice. The boundary conditions here are that of a crystal corner: 
all dimers are aligned the same way deep enough in each of the 3 sectors. On the right is 
(the boundary of) the amoeba of a straight line. 

Most random surface models cannot be solved exactly, and authors are content to prove 
qualitative results about the surface tension, the existence of facets, the set of gradient Gibbs 
measure, etc. 

We will prove in this paper, however, that models based on perfect matchings (on any 
weighted doubly-periodic bipartite graph G in the plane) are exactly solvable in a rather 
strong sense. Not only can we derive explicit formulas for the surface tension — we also 
explicitly classify the set of Gibbs measures on tilings and explicitly compute the local 
probabilities in each of them. These results are a generalization of [2] where similar results 
for G = Z 2 with constant edge weights were obtained. 



3 



The theory has some surprising connections to algebraic geometry. In particular, the 
phase diagram of the Gibbs measures for dimers on a particular graph is represented by 
the amoeba of an associated plane algebraic curve, the spectral curve, see Theorem 14.11 We 
recall that by definition the amoeba of an affine algebraic variety X G C n (plane 
curve, in our case) is the image of X under the map taking coordinates to the logarithms of 
their absolute value. See Figure El for an illustration. The so-called Ronkin function of the 
spectral curve turns out to be the Legendre dual of the surface tension (Theorem 13. 6J) . 

Crystal facets in the model are in bijection with the components of the complement of the 
amoeba. In particular, the bounded ones correspond to compact holes in the amoeba; the 
number of bounded facets equals the genus of the spectral curve. By the Wulff construction, 
the Ronkin function describes the limit height function in the dimer model for suitable 
boundary conditions (which can be interpreted as the shape of a partially dissolved crystal 
corner). For example, the limit shape in the situation shown in Figure ^ is the Ronkin 
function of the straight line. It has genus zero and, hence, has no bounded facets. A more 
complicated limit shape, in which a bounded facet develops, can be seen in Figures 121 and |5] 




-3-2-10 1 2 3 



Figure 2: On the left are the level sets of perimeter 10 or longer of the height function of 
a random volume-constrained dimer conhguration on Z 2 (with 2x2 fundamental domain). 
The height function is essentially constant in the middle — a facet is developing there. The 
intermediate region, in which the height function is not approximately linear, converges to 
the amoeba of the spectral curve, which can be seen on the right. The spectral curve in this 
case is a genus 1 curve with the equation z + z~ l + w + w~ l = 6.25. 

Crystals that appear in nature typically have a small number of facets — the slopes of 
which are rational with respect to the underlying crystal lattice. But laboratory conditions 
have produced equilibrium surfaces with up to sixty identifiably different facet slopes [T7j . 
It is therefore of interest to have a model in which it is possible to generate crystal surfaces 
with arbitrarily many facets and to observe precisely how the facets evolve when weights 
and temperature are changed. 
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For another surprising connection between dimers and algebraic geometry see [TEj . 
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2 Definitions 

2.1 Combinatorics of Dimers 

2.1.1 Periodic bipartite graphs and matchings 

Let G be a Z 2 -periodic bipartite planar graph. By this we mean G is embedded in the plane 
so that translations in Z 2 act by color-preserving isomorphisms of G - isomorphisms which 
map black vertices to black vertices and white to white. An example of such graph is the 
square-octagon graph, the fundamental domain of which is shown in Figure El More familiar 




Figure 3: The fundamental domain of the square-octagon graph. 

(and, in a certain precise sense, universal) examples are the standard square and honeycomb 
lattices. Let G n be the quotient of G by the action of nl?. It is a finite bipartite graph on 
a torus. 

Let M. (G) denote the set of perfect matchings of G. A well-known necessary and sufficient 
condition for the existence of a perfect matching of G is the existence of a unit flow from 
white vertices to black vertices, that is a flow with source 1 at each white vertex and sink 
1 at every black vertex. If a unit flow on G exists, then by taking averages of the flow over 
larger and larger balls and subsequential limits one obtains a unit flow on G±. Conversely, 
if a unit flow on G\ exists, it can be extended periodically to G. Hence G\ has a perfect 
matching if and only if G has a perfect matching. 

2.1.2 Height function 

Any matching M of G defines a white-to-black unit flow u: flow by one along each matched 
edge. Let M be a fixed periodic matching of G and ujq the corresponding flow. For any 
other matching M with flow lu, the difference u — ujq is a divergence-free flow. Given two 
faces fo, fi let 7 be a path in the dual graph G* from f to f\. The total flux of u> — u across 
7 is independent of 7 and therefore is a function of f\ called the height function of M. 
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The height function of a matching M is well-defined up to the choice of a base face fo and 
the choice of reference matching M . The difference of the height functions of two matchings 
is well-defined independently of M . 

A matching Mi of G\ defines a periodic matching M of G; we say Mi has height change 
(j, k) if the horizontal and vertical height changes of M for one period are j and k respectively, 
that is 

h(v + (x, y)) = h{v) + jx + ky 

where h is the height function on M. The height change is an element of Z 2 , and can be 
identified with the homology class in Hi(T 2 , Z) of the flow uj\ — ujq. 

2.2 Gibbs measures 

2.2.1 Definitions 

Let £ be a real- valued function on edges of G±, the energy of an edge. It defines a periodic 
energy function on edges of G. We define the energy of a finite set M of edges by S(M) = 

A Gibbs measure on M.(G) is a probability measure with the following property. If we 
fix the matching in an annular region, the matchings inside and outside the annulus are 
independent of each other and the probability of any interior matching M is proportional to 
e~ £ ( M \ An ergodic Gibbs measure (EGM) is a Gibbs measure on JA(G) which is invariant 
and ergodic under the action of Z 2 . 

For an EGM /x let s = K[h(v+(l,0))-h(v)] and t = E[h(v+(0,l))-h(v)) be the expected 
horizontal and vertical height change. We then have E[/i(f + (x,y)) — h(v)} = sx + ty. We 
call (s, t) the slope of /i. 

2.2.2 Gibbs measures of fixed slope 

On M.{G n ) we define a probability measure /j n satisfying 

e -S{M) 
Hn{M) = , 

for any matching M G M.(G n ). Here Z is a normalizing constant known as the partition 
function. 

For a fixed (s,t) G M 2 , let JA s ,t{G n ) be the set of matchings of G n whose height change 
is (L^sJ, L n ^J)- Assuming that A4 S;t (G n ) is nonempty, let /i n (s,t) denote the conditional 
measure induced by fi n on M. s t (G n ). 

The following results are found in Chapters 8 and 9 of [TH] : 

Theorem 2.1 ( [18J ) . For each (s,t) for which Ai s j{G n ) is nonempty for n sufficiently 
large, fi n (s,t) converges as n — > oo to an EGM fi(s,t) of slope (s,t). Furthermore fi n itself 
converges to /i(so,to) where (so,t ) is the limit of the slopes of //„. Finally, if {so,t ) lies in 
the interior of the set of (s,t) for which A4 s ,t{G n ) is nonempty for n sufficiently large, then 
every EGM of slope (s, t) is of the form //(s, t) for some (s, t) as above; that is, /i(s, t) is the 
unique EGM of slope (s,t). 
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2.2.3 Surface tension 

Let 

Z s , t (Gn) = £ e- £ ^ 
MeM a ,t{Gn) 

be the partition function of M. s t (G n ). Define 

Z Sit (G) = lim Z Sjt (G n ) 1/n2 . 

n— >oo 

The existence of this limit is easily proved using subadditivity as in |2j. The function Z Sjt (G) 
is the partition function per fundamental domain of /x(s, t) and 

<7(s,f) = -logZ a>t (G) 

is called the surface tension or /ree energy per fundamental domain. The explicit form of 
this function is obtained in Theorem 13.61 

The measure /x(so^o) i n Theorem 12.11 above is the one which has minimal free energy 
per fundamental domain. Since the surface tension is strictly convex (see Chapter 8 of ^H] 
or Theorem 13 . 71 below) . the surface-tension minimizing slope is unique and equal to (s ,t ). 

2.3 Gauge equivalence and magnetic field 
2.3.1 Gauge transformations 

Since G is bipartite, each edge e = (w, b) has a natural orientation: from its white vertex 
w to its black vertex b. Any function / on the edges can therefore be canonically identified 
with a 1-form, that is, a function on oriented edges satisfying /(— e) = —/(e), where — e 
is the edge e with its opposite orientation. We will denote by Q}(Gi) the linear space of 
1-forms on G\. Similarly, Q° and Q 2 will denote functions on vertices and oriented faces, 
respectively. 

The standard differentials 

o a> n 1 ± n 2 - o 

have the following concrete meaning in the dimer problem. Given two energy functions Si 
and £ 2) we sa Y th & t they are gauge equivalent if 

£i = £ 2 + df, fen , 

which means that for every edge e = (w, b) 

£ 1 (e)=£i(e)+/(b)-/(w) l 

where / is some function on the vertices. It is clear that for any perfect matching M, the 
difference Si(M) — £z(M) is a constant independent of M, hence the energies £-\ and £2 
induce the same probability distributions on dimer configurations. 
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2.3.2 Rotations along cycles 

Given an oriented cycle 

7 = {w , b , wi, bi, . . . , b fc _i, w k } , w k = \N , 
in the graph G\, we define 

k-l 

£ = [^K- b i) - £( w *+i» b i)] • 

i=l 

It is clear that £x and £ 2 are gauge equivalent if and only if f S\ = £2 for all cycles 7. We 
call f £ the magnetic flux through 7. It measures the change in energy under the following 
basic transformation of dimer configurations. 

Suppose that a dimer configuration M is such that every other edge of a cycle 7 is 
included in M. Then we can form a new configuration M 1 by 

M' = M A 7 , 

where A denotes the symmetric difference. This operation is called rotation along 7. It is 
clear that 

£(M') =£{M)± j £. 

The union of any two perfect matchings M\ and M 2 is a collection of closed loops and one 
can obtain M 2 from Mi by rotating along all these loops. Therefore, the magnetic fluxes 
uniquely determine the relative weights of all dimer configurations. 

2.3.3 Magnetic field coordinates 

Since the graph G\ is embedded in the torus, the gauge equivalence classes of energies are 
parametrized by R F_1 ©IR 2 , where F is the number faces of G\. The first summand is d£, a 
function on the faces subject to one relation: the sum is zero. We will denote the function 
d£ G Vt 2 (Gi) by B z . The other two parameters 

(B x , B y ) E R 2 

are the magnetic flux along a cycle winding once horizontally (resp. vertically) around the 
torus. 

In practice we will fix B z and vary B x ,B y , as follows. Let 7^ be a path in the dual of 
G\ winding once horizontally around the torus. Suppose that k edges of G\ are crossed by 
7 X . On each edge of G\ crossed by j x , add energy ±^AB X according to whether the upper 
vertex (the one to the left of 7^ when 7^ is oriented in the positive x-direction) is black or 
white. Similarly, let 7^ be a vertical path in the dual of Gi, crossing k' edges of Gi, add 
±jjAB y to the energy of edges crossed by 7 y according to whether the left vertex is black 
or white. 

The new magnetic field is now B' = B + (0, AB X , AB y ). This implies that the change 
in energy of a matching under this change in magnetic field depends linearly on the height 
change of the matching: 
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Lemma 2.2. For a matching M of Gi with height change (h x , h y ) we have 
£ B ,{M) - e B >(M Q ) = S B (M) - £ B (M ) + AB x h x + AB y h y . 

3 Surface tension 

3.1 Kasteleyn matrix and characteristic polynomial 

3.1.1 Kasteleyn weighting 

A Kasteleyn matrix for a finite bipartite planar graph T is a weighted, signed adjacency 
matrix for T, whose determinant is the partition function for matchings on T. It can be 
defined as follows. Multiply the edge weight of each edge of T by 1 or —1 in such a way that 
the following holds: around each face there are an odd number of — signs if the face has 
mod 4 edges, and an even number if the face has 2 mod 4 edges. This is always possible 
0. 

Let K = (i^ w b) be the matrix with rows indexing the white vertices and columns indexing 
the black vertices, with K w b being the above signed edge weight ±e _£ ^ w ' b ^ (and if there 
is no edge). Kasteleyn proved [Hj that | deti^j is the partition function, 

\detK\ = ZiY) = e ~ £(m) - 
meM(r) 

3.1.2 Periodic boundary conditions 

For bipartite graphs on a torus T 2 , one can construct a Kasteleyn matrix K as above 8J. 
Then | det K | is a signed sum of weights of matchings, where the sign of a matching de- 
pends on the parity of its horizontal and vertical height change. This sign is a function on 
H\(T 2 , Z/2), that is, matchings with the same horizontal and vertical height change modulo 
2 appear in det K with the same sign. Moreover of the four possibly parity classes, three 
have the same sign in det K and one has the opposite sign [19 . In other words, this sign 
function is one of the 4 spin structures or theta characteristics on the torus. 

Which one of the 4 spin structures it is depends on the choices in the definition of the 
Kasteleyn matrix. By an appropriate choice we can make the (0, 0) parity class (whose 
height changes are both even) have even sign and the remaining classes have odd sign, that 
is, det K = Mqo — Miq — Mqi — Mn, where Moo is the partition function for matchings with 
even horizontal and vertical height changes, and so on. 

The actual partition function can then be obtained as a sum of four determinants 

z = + + Z^ + Z^), 

where Z^ is the determinant of K in which the signs along a horizontal dual cycle (edges 
crossing a horizontal path in the dual) have been multiplied by (— I) 6 and along a vertical 
cycle have been multiplied by (— 1) T . (Changing the signs along a horizontal dual cycle 
has the effect of negating the weight of matchings with odd horizontal height change, and 
similarly for vertical.) For details see |51 119j. 
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3.1.3 Characteristic polynomial 

Let K be a Kasteleyn matrix for the graph G\ as above. Given any positive parameters z 
and w, we construct a "magnetically altered" Kastelyn matrix K(z, w) from K as follows. 

Let 7^ and j y be the paths introduced in Section 12.3.31 Multiply each edge crossed by 
'-fx by z ±l depending on whether the black vertex is on the left or on the right, and similarly 
for 7 r See Figure 0]for an illustration of this procedure in the case of the honeycomb graph 
with 3x3 fundamental domain. We will refer to P(z, w) = det K(z, w) as the characteristic 
polynomial of G. 




Figure 4: The operator K(z, w) 



For example, for the square-octagon graph from Figure |H] this gives 

P(z, w) = z H \- w -\ h 5 . 

z w 

Recall that Mq denotes the reference matching in the definition of the height function and 
ujq denotes the corresponding flow. Let x denote the total flow of ujq across -y x and similarly 
let 2/0 the total flow of u) Q across j y . The above remarks imply the following: 

Proposition 3.1. We have 

P(z,w) = z- xo w- yo e- £{M) z h *w hy (-l) hxh * +h * +h » 
MeM(Gi) 

where h x = h x (M) and h y = h y (M) are the (integer) horizontal and vertical height change 
of the matching M and £(M) is its energy. 

Since G\ has a finite number of matchings, P(z, w) is a Laurent polynomial in z and 
w with real coefficients. The coefficients are negative or zero except when h x and h y are 
both even. Note that if the coefficients in the definition of P were replaced with their 
absolute values (i.e., if we ignored the (— \^ h ^ h y +hx+h v factor), then P(l, 1) would be simply 
the partition function Z(G\), and P(z,w) (with z and w positive) would be the partition 
function obtained using the modified energy S'(M) = S(M) + h x log 2 + h y \ogw. With 
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the signs, however, P((— l) e , (— l) r ) = and the partition function may be expressed in 
terms of the characteristic polynomial as follows: 

Z = \ (~P(1, 1) + P(l, -i) + P(~l, 1) + p(-i, -i)) 

As we will see, all large-scale properties of the dimer model depends only on the polyno- 
mial P(z, w). 

3.1.4 Newton polygon and allowed slopes 

By definition, the Newton polygon N(P) of P is the convex hull in IR 2 of the set of integer 
exponents of monomials in P, that is 

N(P) = convex hull {(j, k) G Z 2 | z^w k is a monomial in P{z, w)} . 

Proposition 3.2. The Newton polygon is the set of possible slopes of EGMs, that is, there 
exists an EGM fj,(s, t) if and only if (s, t) G N(P). 

Proof. First we prove that (s, t) must lie in N(P). A translation-invariant measure of average 
slope (s, t) determines a unit white-to-black flow on G\ with vertical flux s and horizontal 
flux t: the flow along an edge is the probability of that edge occurring. However the set 
of matchings of G 1 are the vertices of the polytope of unit white-to-black flows of G 1: and 
the height change (s, t) is a linear function on this polytope. Therefore (s, t) is contained in 
N(P). 

The absolute value of the coefficient of z l w^ in P is the weighted sum of matchings of 
G\ with height change Therefore there is a matching corresponding to each slope at 

the corners of N(P). Next, we show that given two periodic matchings of slopes (si,ti) and 
(s 2 ,t 2 ), there is a periodic matching (possibly with larger period) with their average slope. 
Suppose the two matchings have the same period (otherwise take the 1cm of their periods). 
The union of the two matchings is a set of cycles; one can change from one matching to 
the other by rotating along each cycle. If the slopes (si,£i) and (s2,t 2 ) are unequal some 
of these cycles necessarily have nonzero homology the torus, so that rotating along them 
will change the slope. On a fundamental domain of twice the size, shift one of each pair of 
all such cycles; this way one creates a matching with slope which is the average of the two 
initial slopes. Taking limits of these periodic constructions, one can produce a tiling with 
asymptotic slope corresponding to any w G N(P). □ 

Note that changing the reference matching Mq in the definition of the height function 
merely translates the Newton polygon. 

3.2 Asymptotics 

3.2.1 Enlarging the fundamental domain 

Characteristic polynomials of larger graphs may be computed recursively as follows: 
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Theorem 3.3. Let P n be the characteristic polynomial of G n . Then 

Zq=Z, Wq=W 

Proof. We follow the argument of |2j where this fact is proved for grid graphs. Since sym- 
metry implies that the right side is a polynomial in z and w, it is enough to check this 
statement for positive values of z and w. View the Kastelyn matrix K n (z,w) of G n as a 
linear map from the space V w of functions on white vertices of G n to the space V& of functions 
on black vertices. When a and j3 are nth roots of unity, let and V^'^ be the subspaces 
of functions for which translation by one period in the horizontal or vertical direction corre- 
sponds to multiplication by a and (3 respectively. Clearly, these subspaces give orthogonal 
decompositions of V w and H, and K n (z,w) is block diagonal in the sense that it sends an 
element in V^ ,/3 to an element in Vf? . We may thus write det K n (z, w) as a product of 
the determinants of the n 2 restricted linear maps from to V^ Q,/3 ; these determinants are 
given by det K(az l / n } I3w l ' n ). □ 

This recurrence relation allows us to compute partition functions on general G n in terms 
of P: 

Corollary 3.4. 

z(G n ) = 1 -(-zr+zr ) +z^+z^), a) 

where 

zt ) =Pn((-i) 9 ,(-m= n n p M- ( 2 ) 

z n = (-l) e , w n = (-l) T 



3.2.2 Partition function per fundamental domain 

We are interested in the asymptotics of Z(G n ) when n is large. The logarithm of the 
expression (J2J is a Riemann sum for an integral over the unit torus T 2 = {(z, w) G C 2 : \z\ — 
\w\ = 1} of logP; the evaluation of this is a Riemann sum so we have 

1 , 1 f , ^/ .dzdw 

n^ Z '« , = wI l0EP( ^ , ~ + °( 1 ) 

on condition that none of the points 

{(z,w):z n =(-l) e ,w n = (-iy} (3) 

falls close to a zero of P. If it does, such a point will affect the sum only if it falls within 
e -o(n ) Q f a zero Q f p ^ an( j j n an y case can on iy decrease the sum). In this case, for any 

n' near but not equal to n, no point of the form Q will fall so close to this zero of P. In 
Theorem 15 . 1 1 below we prove that P has at most two simple zeros on the unit torus. It follows 
that only for a very rare set of n does the Riemann sum not approximate the integral. 
Since the partition function Z(G n ) satisfies 

Zt)< Z (G n )<2m a x{Z^}, 
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we have 

lim'ilogZ(Gy = / \ogP(z,w)-—, 

where the lim' means that the limit holds except possibly for a rare set of ns. But now a 
standard subadditivity argument (see e.g. 0) shows that Z(G n ) 1 ^ n < Z{G m ) l ^ m (l + o(l)) 
for all large m so that in fact the limit exists without having to take a subsequence. 

Theorem 3.5. Under the assumption that P(z,w) has a finite number of zeros on the unit 
torus T 2 , we have 

logZ^ lim - 2 \ogZ(G n ) = -i- f logP(z,w)- — - 

The quantity Z is the partition function per fundamental domain. 

3.2.3 The amoeba and Ronkin function of a polynomial 

Given a polynomial P(z, w), its Ronkin function is by definition the following integral 

F ^y) = 7^l XogP^z^w)-—. (4) 

A closely related object is the amoeba of the polynomial P defined as the image of the curve 
P(z, w) = in C 2 under the map 

(z, w) i — > (log \z\, log |w|) . 

We will call the curve P(z,w) = the spectral curve and denote its amoeba by A(P). 

It is clear that the integral (@J) is singular if and only if (x, y) lies in the amoeba. In fact, 
the Ronkin function is linear on each component of the amoeba complement and strictly 
convex over the interior of the amoeba (in particular, implying that each component of 
M 2 \ A(P) is convex). This and many other useful facts about the amoebas and Ronkin 
function can be found in |14j . See Figures 03 and |H1 for an illustration of these notions. 

We distinguish between the unbounded complementary components and the bounded 
complementary components. 

3.2.4 Surface tension 

Theorem 12.11 gave, for a fixed magnetic field, a two-parameter family of EGMs {/i(s,t)}. 
Let us vary the magnetic field as in Section 12.3.31 Let fi n (B x ,B y ) be the measure on dimer 
configurations on G n in the presence of an additional parallel magnetic field B x , B y . We 
define fi(B x , B y ) to be the limit of fj, n (B x , B y ) as n — > oo and let Zb x ,b v be its partition 
function per fundamental domain. 

We can compute Zb x ,b v in two different ways. On the one hand, the characteristic 
polynomial becomes P(e Bx z, e By w) and, hence, by Theorem l3.3l we have Zb x ,b v = F(B X , B y ), 
where F is the the Ronkin function of P. On the other hand, using Lemma 12.21 and basic 
properties of the surface tension (see Section 12. 2. 3j) we obtain 

F(B X , B y ) = max (-a(s, t) + sB x + tB y ) . (5) 

(s.t) 
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Figure 5 : The curved part of minus the Ronkin function ofz + ^ + w+^ + 5. This is the limit 
height function shape for square-octagon dimers with crystal corner boundary conditions. 




Figure 6: The amoeba corresponding to the model of Figure^ Its complement has one 
bounded component for each of the five interior integer points of N(P) = {x : \x\ < 2}. 
It also has two "semi-bounded" (i.e., contained in a strip of finite width) components, cor- 
responding to two of the non-corner integer points on the boundary of N(P). The four 
large components correspond to the corner vertices of N(P). In the case of equal weights, 
the holes in the amoeba shrink to points and only four large unbounded components of the 
complement are present. 
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In other words, F is the Legendre dual of the surface tension. Since the surface tension is 
strictly convex, the Legendre transform is involutive and we obtain the following 

Theorem 3.6. The surface tension a(s,t) is the Legendre transform of the Ronkin function 
of the characteristic polynomial P. 

Recall that the Ronkin function is linear on each component of the amoeba complement. 
We will call the corresponding flat pieces of the graph of the Ronkin function facets. They 
correspond to conical singularities (commonly referred to as "cusps" ) of the surface tension 
a. The gradient of the Ronkin function maps M 2 to the Newton polygon N(P). It is known 
that the slopes of the facets form a subset of the integer points inside N(P) ^3]. Therefore, 
we have the following immediate corollary: 

Corollary 3.7. The surface tension a is strictly convex and is smooth on the interior of 
N(P), except at a subset of points in 7? fl N(P). Also, a is a piecewise linear function on 
dN(P), with no slope discontinuities except at a subset of points in Z 2 fl dN(P). 

In Section we will see that the spectral curves of dimer models are always very spe- 
cial real plane curves. As a result, their amoebas and Ronkin function have a number of 
additional nice properties, many of which admit a concrete probabilistic interpretation. 

The following Figure shows the Legendre dual of the Ronkin function from Figure 
It is the surface tension function for certain periodically weighted dimers on the square 
grid with 2x2 fundamental domain and also for the uniformly weighted dimers on the 
square-octagon graph. 




Figure 7: (negative of) Surface tension for the square-octagon graph 



4 Phases of the dimer model 
4.1 Frozen, liquid, and gaseous phases 

We will show that EGMs with distinctly different qualitative properties are possible in a 
general periodic dimer model. The different types of behavior can be classified as frozen, 
liquid, and gaseous. We will take the fluctuation of the height function as the basis for the 
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classification. This, as it will turn out, is equivalent to the classification by the rate of decay 
of correlations. 

Let / and /' be two faces of the graph G and consider the height function difference h(f) — 
h(f'). An EGM is called a frozen phase if some of the height differences are deterministic — 
i.e., there exist distinct / and /' arbitrarily far apart for which h(f) — h(f') is deterministic. 
An example is the delta-measure on the brick-wall matching of the square grid. 

A non-frozen EGM /i is called a gaseous phase or smooth phase if the height fluctuations 
have bounded variance, i.e., the /i variance of the random variable h(f) — h(f') is bounded 
independently of / and /'. A non-frozen EGM /i is called a liquid phase or rough phase if 
the /i-variance of the height difference is not bounded. The difference between the smooth 
and rough phases is illustrated in Figure |H1 




Figure 8: All cycles of length ten or longer in the union of two random perfect matching of 
Z 2 with 4x4 fundamental domain. The weight of one edge equals to 10 and all other edges 
have weight 1. The amoeba for this case is plotted in Figure® The slope on the left is (0, 0) 
and this is a smooth phase. The rough phase on the right has slope (0, 0.5). 

We will prove in Theorem 14.51 that in the liquid phase the variance of h(f) — h(f') grows 
universally like 7r _1 times the logarithm of the distance between / and /'. The following is 
our main result about phases: 

Theorem 4.1. The measure n(B x ,B y ) is frozen, liquid, or gaseous, respectively, when 
(B x ,B y ) is respectively in the closure of an unbounded complementary component of A(P), 
in the interior of A(P), or in the closure of a bounded complementary component of A(P). 

This theorem is proved in the next three sections. In Corollary 9.1.2 of JB] there is a 
different proof that when (s,t) lies in the interior of N(P), /i(s,t) can only be smooth if s 
and t are integers. 
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We will see that in the liquid and gaseous phases the edge-edge correlations decay poly- 
nomially and exponentially, respectively. In the frozen case, some edge-edge correlations do 
not decay at all. 

4.2 Frozen phases 

4.2.1 Matchings and flows 

Recall the interpretation of a matching as a black-to-white unit flow. If M is a matching and 
Mo is the reference matching in the definition of the height function, then the difference of the 
flows M — Mo defines a divergence-free flow. The height function of M is the corresponding 
flux, that is, for two faces fx, f 2 , h(f 2 ) — h(fx) is the amount of flow crossing any dual path 
from fx to / 2 . For two adjacent faces fx, f 2 let d(fi, f 2 ) be the maximal possible (oriented) 
flow along the edge e between them (where e is oriented so that fx is on its left). This 
is the forward capacity of the oriented edge e. That is, if e ^ Mo, its capacity is 1 from 
its black vertex to its white vertex, and in the other direction; if e G Mo, its capacity 
is 1 from its white vertex to its black vertex, and zero in the reverse direction. For any 
two faces fx and f 2 let D(fx, f 2 ) be the minimum, over all dual paths from fx to f 2 , of the 
sum of the capacities of the segments oriented to cross the path from left to right. By the 
max-flow-min-cut theorem, a function h is the height function for a tiling if and only if 

V/i,/ 2 D{fx,h)>h{f 2 )-h{fx). 

See j2Ql Oj for a reference. 

Now let (s,t) G R 2 . If there is no tiling with height function having slope (s, t) then there 
is a face / and (x, y) G Z 2 such that D(f, f+(x, y)) < sx+ty. We claim that in this case there 
is a face path from / to some translate / + (x', y') on which D(f,f+ (x', y')) < sx' + ty' and 
all faces along this path are of distinct types, that is, are not translates of each other (except 
for the first and last faces). To see this, note that if a face path fx, f 2 , ■ ■ ■ , fk passes through 
two faces of the same type, say fi and fj, then one of the two paths fx, ■ ■ ■ , fi, fj+i, ■ ■ ■ , fk 
and fi, . . . , fj will necessarily satisfy the strict inequality. 

But up to translation there are only a finite number of face paths which start and end 
at the same face type and which pass through each face type at most once. Each such path 
gives one restriction on the slope: D(f\, f 2 ) > sx + ty where (x, y) is the translation of the 
path. 

In particular the Newton polygon N(P) is the set of (s, t) defined by the intersection 
of the inequalities {(s, t) \ sx + ty < D(fx, ^2)}, one for each of the above finite number of 
paths. If (s,t) is on the edge of N(P), the path 7 from the corresponding inequality has 
maximal flow, that is, in a tiling of slope (s, t) all edges on 7 are determined: they must 
occur with probability 1 or 0. 

4.2.2 Frozen paths 

When (B x , B y ) is in an unbounded component of the complement of the amoeba, we prove 
that n(B x ,B y ) is in a frozen phase. 

The slope (s, t) of [i(B x , B y ) is an integer point on the boundary of N(P). By the argu- 
ment of the previous section, there is a face path 7 on Gx, with homology class perpendicular 
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to (s,t), for which every edge crossing each lift of 7 is present with probability 1 or for 
fi(B x ,B y ). These lifts constitute frozen paths in the dual G'. Edges which are in different 
components of the complement of the set of frozen paths are independent. 

For each corner of N(P) there are two sets of frozen paths, with different asymptotic 
directions. The components of the complements of these paths are finite sets of edges. The 
edges in each set are independent of all their translates. 

4.3 Edge-edge correlations 

For a finite planar graph T the inverse of the Kasteleyn matrix determines the edge prob- 
abilities: the probability of a set of edges {ei, . . . , e^} being in a random matching is the 
determinant of the corresponding submatrix of K' 1 , times the product of the edge weights 



For a graph on a torus the corresponding statement is more complicated: We have 

Theorem 4.2 (|2J). The probability of edges {e\ = (w 1; t>i), . . . , ej. — (w^, b^)} occurring is 
in a random matching of G n equals Y\K(vjj, ^j) times 



where Q(z, w) is one of a finite number of polynomials (depending on where w and b sit in 
their respective fundamental domains: Q/P is an entry of K{z,w)~ 1 where K(z,w) is the 
matrix of Lemma I3.1J1 and (x, y) G Z 2 is the translation taking the fundamental domain 
containing w to the fundamental domain containing b. 
This expression is a Riemann sum for the integral 



except near the zeros of P. However the contribution for the root (z,w) nearest to a zero 
of P is negligible unless (z,w) is at distance O(-tt) of the zero. But if this is the case then 
replacing n with any n' at distance at least 0(^/n) from n makes the contribution for this 
root negligible. Thus we have 



ma. 




Here the determinants det(ii"^ r 1 (bj, w^)) are k x k minors of KZ. . The asymptotics of 
this expression are again complicated by the zeros of ? on T 2 . The entries in have the 
form (see 0) 





lim'^ r (w,b) 



1 f Q(z,w)w x z y dw dz 



(27ii) 2 Jj2 P(z, w) w z 
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where the limit is taken along a subsequence of n's. 

Since all the K~\ T have the same limit along a subsequence of ns, their weighted average 
(as in ©) with weights ±Zg T /2Z (weights which sum to one and are bounded between — 1 
and 1) has the same (subsequential) limit. This subsequential limit defines a Gibbs measure 
on A4(G). By Theorem 12. 1| this measure is the unique limit of the Boltzmann measures on 
G n . Thus we have proved 

Theorem 4.3. For the limiting Gibbs measure fi = lim^oo/in, the probability of edges 
{ei, . . . , ei} where tj = (wj, bf), is 

^JJ K(wj, bj^j det(K~\b k , Wj))i< jt k<i , 

where, assuming b and w are in a single fundamental domain, 

K-\b,w+(x,y)) = -±- [ K-\z,w) bvj w x zy — -. (7) 
We reiterate that K~ 1 (z,w)bw = Qbw(z, w)/P(z, w), where Qbw is a polynomial in z and 

w. 

4.4 Liquid phases (rough non-frozen phases) 
4.4.1 Generic case 

When (B x , B y ) is in the interior of the amoeba, P(e Bx z, e By w) either has two simple zeros on 
the unit torus or a real node on the unit torus (see Theorem l5.11 below). In the case of simple 
zeros (see Lemma below), K~ x (b, w) decays linearly but not faster, as \w — b\ — > oo. This 
implies that the edge covariances decay quadratically: 

Cov(ei, e 2 ) := Pr(d and e 2 ) - Pr(ei) Pr(e 2 ) 

= -K( Wl , b 1 )K(w 2) b 2 )K-\b 2 , yN 1 )K-\b 1 , w 2 ) . 

In section 14.4.21 we show that in the case of a real node we have similar behavior. 

Lemma 4.4. Suppose that \z \ = \w \ = 1, Im(— j3wo/az ) > 0, x,y G Z, and R(z,w) is a 
smooth function on T 2 with a single zero at (zq,Wq) and satisfying 

R(z, w) = a(z — z ) + f3{w — w ) + 0(\z — z \ 2 + \w — w \ 2 ). 

Then we have the following asymptotic formula for the Fourier coefficients of R^ 1 : 

1 r w x z y dz dw ~ w o z o q f 1 \ 

(27r«) 2 J T 2 R(z, w) z w 2iri(xazQ — yf3wo) \x 2 + y 2 J 

Note that if R has k simple zeros, 1/R can be written as a sum of k terms, each of which 
is of the above form. 
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Proof. Replacing z with e ta zo and w with e tb wo we have 

a(z — zq) + (3(w — wo) + 0{. . . ) = az ia + (Iwoib + 0{. . . ). 

By adding a smooth function (whose Fourier coefficients decay at least quadratically) to 1/R 
we can replace 1/R with 

1 

az ia + (5wQib 

The integral is therefore 

„.,x y fir ptt Mxb+ya) in x ? y f°° C°° p i{xb+ya) 

(2vr; 
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rn i(xb+ya) x V poo poo i(xb+ya) 

/ dadb + 0(...) = p^r / — dadb + 0(...). 

J-Tr pw b + az a (2n) 2 i J_ OD J^^ pw b + az a 



We first integrate over the variable a: the integrand is a meromorphic function of a with a 
simple pole in the upper half plane if Im(— f3w b/az ) > 0. Change the path of integration 
to a path from —N to iV followed by the upper half of a semicircle centered at the origin of 
radius N. The residue theorem then yields 

W Z f°° e i(x-y(3wo/azo)b db = W Z ~ j ~ W Z 



2ho>Zq J 2ixq>Zq i(x — y[3w /az ) 2m(xaz — y[3w ) 

which gives the result. □ 

We now compute the variance in the height function. 

Theorem 4.5. Suppose that the zeros of P on T 2 are simple zeros at (zo,wo) and (zq,wo). 
Let a, (3 be the derivatives of P(z,w) with respect to z and w at (zo,wo). Then the height 
variance between two faces f\ and fi is 

Var[h{h) - h(f 2 )\ = I log -0(/ 2 )| +o(log 10(A) -0(/ 2 )|), 

7T 

where is the linear mapping <j>(x + iy) = xazo — y/3u>o- 

Since is a non-degenerate linear mapping, the above expression for the variance is 
equivalent to ^ log \ fi — f 2 \ +o(log ]/i — f 2 \). However it appears that a slightly finer analysis 
would improve the little-o error in the statement to o(l), so we chose to leave the expression 
in the given form. 

Proof. Define h = h — E(/i). Let fi, f 2 , h,f± be four faces, all of which are far apart 
from each other. We shall approximate (h{fi) — h(f 2 ))(h(f 3 ) — h(f^)). To simplify the 
computation we assume that f\ and f 2 are translates of each other, as well as / 3 and / 4 . Let 
fi — 9i: 92, ■ ■ ■ , 9k — f2 be a path of translates of f\ from f\ to f 2 , with g p+ \ — g p being a 
single step in Z 2 . Similarly let fz — g[, g 2 , ■ ■ ■ , g\ — fi be a path from f 3 to f±. We assume 
that these paths are far apart from each other. 
Then 

fe-i i-\ 

(Hfi) - hfoMHh) - Kh)) ^E^M - - hisQ). (8) 

p=l q=l 
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We consider one element of this sum at a time. There are three cases to consider: when 
9p+i ~ 9p an d g' q +i — g' q are both horizontal, both vertical, and one vertical, one horizontal. 
From Lemma f4.4| for x, y large we have 

r Q(z,w)w x z y dzdw _ 2R / w%z%Q(z ,w Q ) \ | 1 
J P(z,w) z w \2iri(xaz — y(3w ) J x 2 + y 2 

Recall that the matrix Q(z,w) satisfies Q(z,w)K(z,w) = P(z,w) ■ Id. Since (zo,Wq) is 
a simple zero of P, K{zq,wq) has co-rank 1. In particular Q{zq ) wq) must have rank 1. We 
write Q(z ,w ) = UV l where V t K(z , w ) = = K(z ,w )U. 

Let di = (wj, bj) be the edges crossing a "positive" face path 7 from g p to g p +i, that is, a 
face path with the property that each edge crossed has its white vertex on the left. Similarly 
let a'j = (yj'j, b'j) be the edges crossing a positive face path 7' from g' q to g' q+ i- Then 

E((h(g p+1 )-h(g p ))(h(g' q+1 ) ~W q ))) = E (5>* " " a J)) 

= J^E(oiaJ) -E(oi)E(aj.) 
».j 

= £ K ( w " b *)^K' b^^-^fij., wOliT- 1 ^, wj). 
».j 

Assuming these faces g p ,g' q are far apart, this is equal to 

~ 4^ A ^ Wj ' b *) K K> b i) x 

/ wiziuy; | ggO^g | Q / 1 

I xaz — y/3w xaz — yPw \x 2 + y 2 

f w^z^u^ | grggg / 1 \\ 

yxazo-y/^wo xaz -ypwo \x 2 + y 2 J J 

When we combine the cross terms we get an oscillating factor Zq x w^ v or its conjugate, 
which causes the sum (jHJ) of these terms when we sum over p and q to remain small. So the 
leading term for fixed p, q is 

Now we claim that if / 2 — /1 = (1, 0) then 

^ ^(w i? bj)^ = Zq—(zq, ™o) = ^o" 
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and when f 2 — fi = (0, 1) then 

^2 K (Wii bi)UiVi = Wq — {Zq, Wq) = w /3, 

and similarly for f'4 — f'3. To see this, note first that the function K(wi, bj)UiVj is a function 
on edges which is a closed 1-form, that is, a divergence-free flow. In particular the sum 
^2 i K(\N i ,bi)UiVi is independent of the choice of face path (in the same homology class). 
We can therefore assume that the face paths 7, 7' are equal to either 7^ or j y according to 
whether they are horizontal or vertical. Suppose for example fz — fi = (1, 0) and differentiate 
Q(z, w)K(z, w) = P(z, w) ■ Id with respect to z and evaluate at (z , w ): we get 

Q z (z , W )K(z , W ) + Q(Z , W )K z (z , Wq) = P z (z , Wq) ■ Id. 

Applying U from the right to both sides, using K(z ,w )U = and Q(z ,w ) = UV 1 , and 
then multiplying both sides by Zq, this becomes 

UV t z K z {zQ,WQ)U = z P z (z ,Wq)U. 
However z Q K z (z , w ) = Y. lx K {™, b ) w ® b = S; e7 K {™n bi)w* <g> b i? so we get 

U ^ K ( w i' bi)ViUi = z P z (z , w )U 

and since U ^ the claim follows. The same argument works for 7 y . 

Recall that (x, y) was the translation between the fundamental domains containing g p 
and g' q . In the sum (JHJ), let (xi, yi) G Z, 2 be the position of the fundamental domain of g p and 
( x 2, 2/2) that °f 9q- Let 21 = Xiaz — y\(3wQ and z 2 = x 2 az — y 2 j3w . The sum (JHJ) becomes 
(up to lower order terms) 

l_ Re /^ (/2) /^ (/4) d Zl dz 2 

27r " 6 A(/i) A(/a) (^l-^) 2 

where is the linear map (x, y) 1— > x/?wo — yc^^o- This integral evaluates to 

^ R , ( (<i>(h) - ttfi))(<i>(h) - <Kf*)) \ 

2vr 2 8 ^(0( /4 )-0(/ 2 ))(0( /3 )-0( /l ))y'- 

To approximate the height variance cr(h(f 2 ) — h(fx)), let fs be close to /1 and close to / 2 
(but still far enough apart on the scale of the lattice so that the above approximations hold). 
Then as \f 2 — f\\ — > 00 while I/3 — /i| and I/4 — f 2 \ are remaining bounded, the variance is 

1 log |0(A) - 0(/ 2 )| + (iog |0(A) - 0(/ 2 )|). 

□ 
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4.4.2 Case of a real node 



In this section we show how to modify the above proof in case P has a real node. This 
happens when the two simple zeros (z ,wo), (zo,Wq) merge into a single zero at one of the 
four points (±1,±1), and P z = P w = there. In this case the slope of the corresponding 
EGM is integral but the amoeba does not have a complementary component, The component 
is reduced to a point (and is therefore not "complementary"). 

The canonical example of this behavior is the square grid with uniform weights and a 
2x2 fundamental domain; in this case P = 4 + z + z~ l + w + w~ 1 , and there is a real node 
at (z, w) = (—1, —1). 

Since i^ _1 (b,w) is a continuous function of the edge weights, so is the variance of the 
height between two points. When P(z,w) has a node, at say (z,w) = (1, 1), the polynomial 
P(z,w) = P(e Bx z,e By w) has two simple zeros on T 2 as long as B x ,B y are sufficiently close 
to but not equal to 0. So the height variance of P is the limit of the height variances of P 
as B x , By — ► 0. 

In fact suppose without loss of generality that the node is at (z,w) = (1, 1) and P has 
the expansion 



where o,i,c6l and . . . denotes terms of order at least 3. Then near the node a point on 
P satisfies either z — I — X(w — 1) + 0(w — l) 2 or z — 1 = X(w — 1) + 0(w — w ) 2 where 
A, A, which are necessarily non-real, are the roots of a + bx + cx 2 = 0. 

The proof of Theorem 14.51 is valid for P except for one assertion, where we ignored the 
cross terms in equation JOJ). Indeed, when z,w are each close to 1 the cross terms do not 
oscillate. However we will show that 2~^i-^( w «5 bi)UiVi = o(\a\ + |/3|) as (z,w) tends to the 
node. Along with the complex conjugate equation this proves that the cross terms make no 
contribution. 

The cross terms of (jHJ) give 



where recall that a = P z , (3 = P w are tending to at the node. We first show that 
J2i K(yji, bi)UiVi = 0(\a\ + \/3\), and similarly for its complex conjugate. Recall the equation 
QK = P ■ Id. Differentiating with respect to z we find 



At a point (z ,wo) on P we have Q = UV* and at (zq,Wq) we have Q = UV t . Applying U 
to the right and evaluating in the limit as (zq,Wq) tends to the node, we have 



so that = V t K z U = '^ /i K(\Ni,bi)UiVi at the node. Since U,V can be chosen polynomial 
in z, w the quantity ^-^(wj, b^UiVt necessarily vanishes to order at least one at the node 
(as do a and (3). 



P(z, w) = a(z - l) 2 + b(z - l)(w - 1) + c(w - l) 2 + . . . , 




(10) 



Q Z K + QK Z = P z ■ Id. 



UV l K z U = 
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Let ci be the limit at the node of ^ J2i -^( w i> bj)C/jVi when — g p — (1, 0) and c 2 the 
same limit when g p+ i — g p — (0, 1), so we may write 



lim — K(yji, bAUiVi = c\dx + c 2 dy 



at the node. The cross terms are then 
1 1 



2-7T 2 \x — y/3/a\ 



((cidxx + c 2 dy l )(c 1 dx 2 + c 2 dy 2 ) + (c 1 dx 1 + c 2 dy 1 )(c l dx 2 + c 2 dy 2 )) . 



A short computation now shows that, since a//3 $ R, this is not a closed 1-form in (xi,y±) 
or (x 2 ,y 2 ) unless Ci = c 2 = 0. 

However since the height function differences h(f\) — h(f 2 ) do not depend on the path 
<7i, . . . , <7fc, the cross terms should necessarily be a closed 1-form. So ci = c 2 = 0. We have 
proved 

Theorem 4.6. Suppose P has a real node (zq,wo) = (±1,±1) on the unit torus, and 

P(z, w) = a(z - z ) 2 + b(z — zq) (w — w ) + c(w - w ) 2 + 

Then the height variance between two faces f\ and f 2 is 

Var[h{h) - h(f 2 )\ = I log -0(/ 2 )| +o(log 10(A) -0(/ 2 )|), 

71 

where is the linear mapping 4>(x + iy) = xz — y\wo, A being the root of a + b\ + cA 2 = 0. 



4.5 Gaseous phases (smooth non-frozen phases) 

When (B x ,B y ) is in a bounded complementary component, P(e~ Bx z, e~ By w) has no zeros 
on the unit torus. As a consequence K _1 (b, w) decays exponentially fast in |b — w|. 

Proposition 4.7. The height variance a(h(fi) — h(f 2 )) is bounded. 

Proof. The height difference h(f\) — h(f 2 )) can be measured along any path from f\ to 
f 2 . Suppose a dual path from f\ to f 2 is chosen so that each edge crosses an edge of G 
with black vertex on its left. Then the height difference is a constant plus the sum of the 
indicator functions of the edges on the path: h(fi) — h(f 2 ) = a« — a«. Suppose we have 
two such paths 71,72, consisting of edges a« and bj respectively, which are close only near 
their endpoints. The height variance is then the sum of the covariances of a« and bj: 

a{h{h) - h{f 2 )) = ^Pr(a, and bj) - Pr(a,) Pr(^). 

However these covariances are exponentially small except when both aj and bj are near the 
endpoints f\ or f 2 . In particular the above summation is a geometric series, which has sum 
bounded independently of the distance between f\ and f 2 . □ 
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4.6 Loops surrounding the origin 



Smooth phases are characterized by their bounded height variance or exponential decay of 
correlation. Here is another characterization of smooth phases in terms of loops in the union 
of two perfect matchings sampled independently. 

Theorem 4.8. A non-frozen EGM fi is smooth if and only if when two perfect matchings M\ 
and M2 are chosen independently from \x, there are almost surely only finitely many cycles 
in Mi U M2 that surround the origin. Similarly, a non-frozen EGM is rough if and only if 
when two perfect matchings Mi and M 2 are chosen independently from \i, there are almost 
surely infinitely many cycles in Mi U M 2 that surround the origin. 

Proof. Since having infinitely many cycles surround the origin is a translation-invariant 
event, it is clear that if \x is an EGM, then there are either /x-almost-surely infinitely many 
cycles or /i-almost-surely finitely many cycles surrounding the origin. If the former is the 
case, it is easy to see that the variance is unbounded; to see this, simply use the fact that, 
conditioned on the positions of the cycles, the two orientations of a given cycle (i.e., which 
alternative set of edges in the cycle belongs to which of the Mj) are equally likely, and orien- 
tations of the cycles are independent of one another. (Note that the orientation of the cycle 
determines whether the height difference of the two height functions goes up or down when 
we cross that cycle.) 

Now suppose that there are almost surely only finitely many cycles surrounding the 
origin. Lemma 8.4.3 of ^H] further implies that if two perfect matchings Mi and M2 are 
sampled independently from /x(s, t), then the union MiUM 2 almost surely contains no infinite 
paths. It follows that the height difference between the two height functions is constant on 
the infinite cluster of faces that are not enclosed in any loops. Lemma 8.3.4 of 18j then 
implies that /x is smooth with respect to the (differently formulated but actually equivalent) 
definition given in Chapter 8 of [18J, and Lemmas 8.1.1 and 8.1.2 imply that the height 
difference variances remain bounded in this case. □ 

5 Maximality of spectral curves 
5.1 Harnack curves 

The characteristic polynomial P(z, w) has real coefficients and, hence, the spectral curve 
P(z, w) = is a real plane curve (in fact, it is more natural to consider the spectral curve as 
embedded in the toric surface corresponding to the Newton polygon of P) . While all smooth 
complex curves of given genus are topologically the same, the number and the configuration 
of the ovals of a real plane curve can be very different. In particular, there is a distinguished 
class of real plane curves, known as Harnack or maximal curves, which have the maximal 
number of ovals (for given Newton polygon) in the, so to speak, best possible position. The 
precise topological definition of a Harnack curve can be found in ^3]; here we will use the 
following alternative characterization of a Harnack curve obtained in ^5]. Namely, a curve 
P(z, w) is Harnack if and only if the map from the curve to its amoeba is 2-to-l over the 
amoeba interior (except for a finite number of real nodes where it is 1-to-l). The main result 
of this section is the following 
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Theorem 5.1. For any choice of nonnegative edge weights the spectral curve P(z,w) = is 
a Harnack curve. 

Harnack curves form a very special and much studied class of curves. Several characteri- 
zations and many beautiful properties of these curves can be found in ^3JE3 (see also [12J). 
We will see that several of them have a direct probabilistic interpretation. 



5.2 Proof of maximality 



Maximality is an important property and several proofs of it are available. In many respects, 
it resembles the notion of total positivity [HE! and the proof given below exploits this analogy. 
Proofs of maximality based on different ideas can be found in ED • 

First observe that by the 2-to-l property, being Harnack is a closed condition, hence it is 
enough to prove that spectral curve is Harnack for a generic choice of weights. Furthermore, 
it is easy to see that we can obtain any periodic planar bipartite graph as a limit case of 
the periodically weighted hexagonal lattice when some of the edge weights are zero. It is 
therefore, enough to consider the case of generic periodic weights on the hexagonal lattice 
with n x n fundamental domain. 




/ 

/ 

















Figure 9: The solution to K(z,w)f = can be constructed layer by layer 



By definition, P(z, w) is a determinant of an n 2 x n 2 matrix K(z, w) whose rows and 
columns are indexed by the n 2 white and black vertices in the fundamental domain. We can 
also write P(z, w) as a determinant of an n x n matrix using transfer matrices as follows. 
The equation det K(z, w) = means that there exists a nonzero function / on black vertices 
annihilated by the operator K(z, w). We can construct such a function row by row as follows: 
given the values of / on a horizontal row of black vertices as in Figure |H1 the equation Kf = 
determines the values of / on the row below it. The corresponding linear map is given by 
—T(w), where T(w) is the transfer matrix of the form 

/ ol\ b\ \ 

a 2 b 2 
a 3 b 3 



T(w) 



\b n w 



h > . 
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Iterating this procedure once around the period, we get a consistency relation, which gives 

P(z,w) = det (z - (-l) n Ti(» • ■■T n (w)) . 

Suppose now that for some point (x, y) G M 2 the torus {\z\ = e x , \ w\ = e y } contains more 
than 2 points of the curve P(z,w). By changing the magnetic field, we can assume that 
y — 0. That is, we can assume that for a pair of points (z\,wi) and (^2,^2) on the spectral 
curve we have 

= \w 2 \ = I , w 2 ^wi,wi, \z 1 \ = \z 2 \. 

Since being not Harnack is an open condition, we can find a nearby curve for which both w± 
and w 2 are roots of unity. It is easy to see (a more general statement is proved in [T2j) that 
we can achieve this by a small perturbation of the dimer weights. So, we can assume that 
w™ = w™ = 1 for some integer m. 

Taking nm as the new horizontal period, we find that the matrix 

M = T 1 (1)---T n (l) (12) 

has more than 2 eigenvalues of the same absolute value. We will now show that this is 
impossible. 

This follows from the following lemma which is a version of a standard argument in the 
theory of total positivity (cf. j3j). 

Lemma 5.2. Suppose that all odd-size minors of a matrix M are nonnegative and that there 
exists k such that all odd-size minors of M k are positive. Then the eigenvalues of M have 
the following form 

Ai > |A 2 | > |A 3 | > |A 4 | > |A 5 | > |A 6 | > |A 7 | > . . . . 

Remark, in particular, that if |A 2 fc| > |A 2 fc + i| then both of these eigenvalues are real. 

Proof. Let us order the eigenvalues of M so that 

I Ai| > |A 2 | > |A 3 | > ... . 

Since all matrix elements of M are nonnegative and all matrix elements of M k are positive, 
the Perron-Frobenius theorem implies that Ai is simple, positive, and that 

Ai > |Aj| , i > 1 . 

Now consider the action of the matrix M in the third exterior power A 3 lR n of the original 
space W n . The matrix elements of this action are the 3x3 minors of M and, hence, Perron- 
Frobenius theorem again applies. The eigenvalues of this action are the numbers 

AjAjAfc , l<i<j<k<n. 

It follows that the number AiA 2 A3 is real, positive, and greater in absolute value than AiA 2 Aj 
for any i > 3. It follows that 

I A3 1 > IA4I . 

Iteration of this argument concludes the proof. □ 

It is immediate to see that any matrix of the form Tj(l) satisfies the hypothesis of the 
lemma. Matrices with nonnegative minors of any given order form a semigroup, so all odd- 
size minors of (fT2*|) are nonnegative. Genericity implies that all odd-size minors of M k are 
positive for some large enough k. This concludes the proof. 
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5.3 Implications of maximality 



5.3.1 Phase diagram of a dimer model 

The 2-to-l property immediately implies that: 

(i) the only singularities of the amoeba map are folds over the boundary of the amoeba; 

(ii) the boundary of the amoeba is the image of the real locus of the spectral curve; 

(iii) to any lattice point in the interior of the Newton polygon corresponds either a bounded 
component of the amoeba complement or an isolated real node of the spectral curve. 
In particular, the number of holes in the amoeba equals the geometric genus of the 
curve. 

For a general plane curve, more complicated singularities of the amoeba map are possible, 
which are then reflected in more complicated singularities of the Ronkin function and, hence, 
of its Legendre dual. In fact, we have used the 2-to-l property in Section|3]in the classification 
of the phases of the dimer model. In this sense, the probabilistic meaning of maximality is 
the absence of any exotic phases with anomalous decay of correlations. 

Part (iii) implies that the gaseous phases persist unless the corresponding component of 
the amoeba complement shrinks to a point and a nodal singularity develops. In particular, 
generically, the spectral curve is smooth and all gaseous phases are present. Conversely, if 
no gaseous phases are present then the spectral curve has the maximal possible number of 
nodes and hence is a curve of genus zero. It is shown in ^2] that the latter case corresponds 
to isoradial dimers studied in jllj . It is also shown in that all Harnack curves arise as 
spectral curves of some dimer model. 

To summarize, part (iii) implies the following: 

Theorem 5.3. The number of gaseous phases of a dimer model equals the genus of the 
spectral curve. For a generic choice of weights, the dimer model has a gaseous phase for 
every lattice point in the interior of the Newton polygon N(P) and a frozen phase for every 
lattice point on the boundary of N(P). 

Part (ii) shows that the phase boundaries can be easily determined. In fact, the amoeba 
of a Harnack curve P(z, w) — can be defined by a single inequality 



where the product is over all 4 choices of signs. Indeed, it is clear that the product in (|T3"|) 
changes sign whenever we cross the amoeba boundary and, generically, is positive when x or 
y are large. We remark that for a general, non-Harnack, curve it is a rather nontrivial task 
to determine its amoeba. 

Using the interpretation of P(±l,±l) as the expectation of (±l) hx (±l) hy (— l) hxh y with 
respect to the measure fii, see Section f3 . 1 . 31 we arrive at the following: 

Theorem 5.4. The minimal free energy measure /i is smooth (that is, frozen or gaseous) if 
and only if the fix-measure of one of the four Hi(T 2 , Z/2) classes of matching s of G\ exceeds 
| ; or can be made more than | by an arbitrarily small perturbation of the dimer weights. 




(13) 
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5.3.2 Universality of height fluctuations 



In Theorem 14. 51 we proved that in a liquid phase the variance of the height function difference 
grows like n- 1 t imes the logarithm of the distance. The proof of that theorem shows that 
the constant in front of the logarithm is directly connected to the number of roots of the 
characteristic polynomial on the unit torus. In particular, maximality was used in the 
essential way to show that this constant is always it^ 1 . In fact, the universality of this 
constant is equivalent to maximality. 

5.3.3 Monge- Ampere equation for surface tension 

It follows from the results of ^H] that the Ronkin function F of a Harnack curve satisfies 
the following Monge- Ampere equation 



for any (x, y) in the interior of the amoeba. By the well-known duality for the Monge- Ampere 
equation, this implies the analogous equation for the surface tension function. 

Theorem 5.5. We have 



It should be pointed out that in [6 J it was argued that for certain class of random surface 
models the equation ([To]) should be satisfied at any cusp of the surface tension. It seems 
remarkable that in our case (|T3jl is satisfies not just at a cusp but identically. For a general 
random surface model, we only expect the left-hand side of (fT5Jl to be positive, by strict 
concavity. 

The geometric meaning of the equations (|14|) and (|15|) is that the gradients of F and a, 
which are mutually inverse maps, are area-preserving, up to a factor. This leads to another 
characterization of Harnack curves as curves with amoebas of maximal possible area for 
given Newton polygon N(P), namely 7r 2 times the area of N(P). It would be interesting to 
find a probabilistic interpretation of this. 

5.3.4 Slopes and arguments 

Recall that by the results of Section 13.2.41 the slope (s,t) of the Ronkin function at a point 
(B x ,B y ) G M 2 equals the slope of the measure fi(B x ,B y ). Suppose that the point (B x ,B y ) 
is in the interior of the amoeba and let (zq, wq) be one of its two preimages in the spectral 
curve. Maximality connects the slope (s, t) with the arguments z and w as follows. 

Theorem 5.6. We have 




(14) 




(15) 



(s,t) = ±— (argu> , argz ) mod Z 2 . 



29 



Proof. We have 



d \ f f , „ „ B B ,dz dw 



1 ' 'if d]ogP y™. 



/p , , log P(e"'z,e^w) 

dB x {2tti) z J J T 2 z w 

2ni J\ w \ =e B y \2m J\ z \ =eBx ° J w 

The inner integral counts, for w fixed, the number of zeros of P(z,w) inside {\z\ = e Bx }. As 
w varies over the unit circle the number of zeros is locally constant with unit jumps whenever 
a zero crosses the circle. These jump points are precisely the points w where P(z, w) has a 
root on T 2 , that is precisely the two points and w where (zo,Wq) and its conjugate are 
the unique zeros of P on T 2 . □ 



6 Random surfaces and crystal facets 

The goal of this section is to review some known facts about random surfaces in order to 
say precisely what Theorem 13 .61 and Theorem 14.11 imply about crystal facets and random 
surfaces based on perfect matchings. We begin with some analytical results. 



6.1 Continuous surface tension minimizers 

A standard problem of variational calculus is the following: given a bounded open domain 
D C M 2 and any strictly convex surface tension function o : R 2 — > R, find the function 
/ : D — ► R. whose (distributional) gradient minimizes the surface tension integral 

1(f) = [ a(Vf(x))dx 

JD 

subject to the boundary condition that / extends continuously to a function / on the 
boundary of D and the volume condition that 

/ f{x) = B 

J D 

for some constant B. In our setting, we may also assume that cr(u) = oo whenever u lies 
outside of the closure of the Newton polygon N(P), so that / is necessarily Lipschitz. The 
following result is well known (see [2], JH] for details and references). 

Proposition 6.1. If there exists any f , satisfying prescribed volume and boundary con- 
straints and satisfying 1(f) < oo, then the surface tension minimizer f is unique and its 
gradient is almost everywhere defined. 

Both Vcr and V/ are functions from R 2 to R 2 . The Euler-Lagrange equation for the 
functional /(/) takes the following form: 
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Proposition 6.2. Let f be the surface tension minimizer described above. Then whenever 
x G D, f is C 2 at x, and a is smooth at f{x), we have: 

div(Vao V/(x)) = C (16) 

for some constant C, which depends on B and /q. If f is also the minimal surface tension 
function when the volume constraint is ignored, then Mb]) holds with C — 0. 

The archetypal solution to the equation in Lemma 16.21 is the Legendre dual of er, which, 
in our setting, is the Ronkin function F. By construction 

Vct o VF = Id 

throughout the amoeba — in particular, the divergence is constant. By analogy with the case 
of the Ronkin function we say that / has a facet of slope u if V/ is equal to u on some open 
subset of D. 

If fo is linear of any slope u on the boundary of D, then it is easy to see that the minimal 
surface tension function, ignoring the volume constraint, is linear of slope u, has (trivially) 
a facet of slope u, and satisfies (fTfiJl with C — 0. If we require, however, that C ^ 0, so that 
the volume constraint exerts some non-zero amount of pressure (upward or downward) on 
the surface, then Lemma 16.21 and Theorem 13.71 imply that the slope of any facet of / must 
be a lattice point inside N(P). In other words, the facet slopes of the Ronkin function F 
represent all possible facet slopes of the dimer model. 

6.2 Concentration inequalities for discrete random surfaces 

In this section, we aim to show that the surface tension minimizing shapes and facets de- 
scribed in Section l6~T1 are approximated by perfect-matching-based discrete random surfaces. 
First, suppose that D is a domain in M 2 , that / is continuous Lipschitz function defined 
on dD, and that / is a surface tension minimizer (with no volume constraint) which agrees 
with f on the boundary. 

Next, denote by -G the infinite graph G whose embedding into M 2 has been re-scaled by 
a factor of 1/n. For example, when G = Z 2 , then -G is a grid mesh that is n times finer than 
G. Suppose that D n is a sequence of simply connected subgraphs of ^G that approximate 
(D, /) from the inside in the sense that 

1. The embedding of D n is contained in D for all n. 

2. The Hausdorff distance between the boundary of D n and the boundary of D tends to 
zero in n. 

3. Each D n admits at least one perfect matching, and for these matchings, the height 
functions h n on the boundary of D n are such that sup|/i„ — /| tends to zero in n 
(where h n is treated as a function on a subset of the points in D, by letting h n (x) 
denote the height at the face of D n containing x). 
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For each n, define v n to be the Boltzmann measure on perfect matchings M n of D n . 
Intuitively, we would expect that for sufficiently large values of n, the normalized function 
h n j n will closely approximate the continuous function / with high u n probability A version of 
this statement is proved in [2] in the case G = Z 2 . Analogous but more general statements (in 
the form of large deviations principles) are discussed in Chapter 7 of fHj (see the paragraphs 
on Lipschitz potentials in Sections 7.3 through 7.5). In both [2] and [TH], the results imply 
that, for a fixed value of e, u n {\h n /n — f \ > e} tends to zero exponentially in n 2 . Also, if 
x G D is a point at which V/(x) = u, we would expect the local statistics of h n , near the 
point x, to be those of the Gibbs measure /i Ul , U2 . More precise versions of this statement can 
be found in Chapter 7 18J . The issue of weighting by enclosed volume is addressed briefly 
in Section 7.5 of [THj . 
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